Edge spin accumulation in semiconductor two-dimensional hole gases 
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The controlled generation of localized spin densities is a key enabler of semiconduc- 
tor spintronics In this work, we study spin Hall effect induced edge spin accumulation 
in a two-dimensional hole gas with strong spin orbit interactions. We argue that it is 
an intrinsic property, in the sense that it is independent of the strength of disorder 
scattering. We show numerically that the spin polarization near the edge induced by 
this mechanism can be large, and that it becomes larger and more strongly localized as 
the spin-orbit coupling strength increases, and is independent of the width of the con- 
ducting strip once this exceeds the elastic scattering mean-free-path. Our experiments 
in two-dimensional hole gas microdevices confirm this remarkable spin Hall effect phe- 
nomenology. Achieving comparable levels of spin polarization by external magnetic 
fields would require laboratory equipment whose physical dimensions and operating 
electrical currents are million times larger than those of our spin Hall effect devices. 

PACS numbers: 75.50.Pp, 85.75.Mm 



INTRODUCTION 



Finite spin densities in semiconductors have tradition- 
ally Q been generated by external magnetic fields, by 
circularly polarized light sources, or by spin injection 
from ferromagnets. Recently there has been consider- 
able interest 2] in an alternate strategy in which edge 
spin densities are generated electrically via the spin Hall 
effect (SHE) 3, 4], i.e., in a planar device by the current 
of spins oriented perpendicular to the plane that is gener- 
ated by and flows perpendicular to an electric fleld. The 
SHE has traditionally been thought of as a consequence 
of spin-dependent chirality in impurity scattering that 
occurs in systems with spin-orbit (SO) coupling 5, 6]. 
Recently it has been recognized that the SHE also has 
an intrinsic contribution due to SO coupling in a perfect 
crystal 7, 8]. The theoretical works ^, 7, 8] and the sub- 
sequent experimental discovery H, 3 of the SHE in both 
strong and weak SO interaction regimes have generated 
substantial interest in both the origins of this unfamiliar 
transport phenomenon and its potential for low power 
consumption electronic devices ]3, 0. The role of dis- 
order in the spin-Hall conductivity is more subtle than 
in more familiar transport coeflicients. Early theoretical 
work reflected a variety of strongly polarized views, all 
the way from a picture with large dissipationless trans- 
verse spin currents T, "s] in which disorder plays no role, 
to demonstrations that the effect is coniple tely eliminated 
by infinitesimally weak disorder |3, [lO: Ull • The theo- 
retical picture is now becoming more clear, if more com- 
plex. It is now established that many of the theoreti- 
cal controversies have arisen because what is seemingly 
the simplest possible model system, a two-dimensional 



electron gas with Rashba SO coupling [T^, turns out to 
have anomalous properties. Although the SHE is absent 
in this model, it is present and robust against disorder 
[13 . [l5l Hg . \u\ in generic systems and in all other models 
that have been speciflcally studied. 

Another challenging issue in SHE theory concerns the 
deflnition of the spin current in SO coupled systems; since 
spin is not conserved the merit of a particular deflnition 
has to be decided on the basis of its utility for the evalua- 
tion of observable properties. The fact that different defl- 
nitions have been proposed which yield spin Hall currents 
that differ even by a sign for some models emphasizes the 
importance of this point 'T^,'T?,'20]. The most important 
observable associated with the SHE is the edge spin accu- 
mulations to which itgives rise; indeed experiments that 
discovered the SHE ^, ^ did so by measuring edge spin 
accumulations. These measurements have partly steered 
the focus of theoretical studies from the experimentally 
ambiguous spin Hall currents to spin Hall accumulations 
21, 22, 23, 24]. The opportunity to directly compare the- 
oretical and experimental values of a SHE related phys- 
ical quantity has been appealing, yet physicists have en- 
tered this territory armed with no analytical theory that 
is able to predict the sign, magnitude, and spatial ex- 
tent of the SHE spin accumulation, at least in the most 
interesting case of strong SO interactions. Particularly 
intriguing in this context is the dual role of SO coupling 
which generates the spin-polarization of carriers while, at 
the same time, causing spin decoherence upon scattering 
off channel edges or off impurities. In the present paper 
we compare a numerical study and experimental results 
of the SHE spin accumulation in two-dimensional hole 
gases (2DHGs). The systems we study are in the diffu- 
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a) spin density profile 



b) spin current profile 



c) edge spin density 




d) bulk spin current 




FIG. 1: Spin Hall effect in a 2DHG. Spatial profile of the 
^-component of spin (a) and of the spin current (b) across 
the channel width (x-direction) . Disorder dependence of the 
accumulated edge spin density (c) and of bulk spin current 

(d). 



sive transport regime in which the sample dimensions are 
large compared to scattering mean-free-path. 



II. THEORY 

First we discuss our theoretical results. The calcu- 
lations are based on the following model Hamiltonian 
jl^,^^ of strongly confined 2DHG at semiconductor het- 
erojunction with structural inversion asymmetry: 



Hr = 



-i\k\ 



(1) 



where m is the carrier effective mass and the momenta 
k = {kl^ A:2)V2 and k± = k^ ± iky. The 2DHG strip 
is defined by hard wall barriers separated by a distance 
Lx = L in the transverse x-direction and by a longitu- 
dinal length Ly = and periodic boundary condition 
in the ^/-direction. Disorder is modeled assuming un- 
correlated, short-range spin-independent scatterers, i.e., 
V{t) = - ^i)' This "white-noise" disorder 

potential satisfies {V{q)V{q'y) = NiV'^6{q_ — q'), where 
V{q) is the Fourier component of V{r) and A^^ = n^L^ is 
the number of impurities; ~ n where n is the 2DHG 
density. The strength of disorder can be characterized by 
the corresponding Born approximation scattering rate, 
h/r = 27rniF^pF, where pp is the 2D density of states 
in the absence of SO coupling. In the following we relate 
all energies to the Fermi energy, Ep, of the 2DHG in the 
absence of SO coupling and disorder, and all lengths to 
the corresponding inverse Fermi wavevector kp^. In our 
studies, Ef is always the largest relevant energy scale. 

The SHE spin-accumulation is generated near the 
edges of the system by a charge current in the ^/-direction. 
Linear response theory predicts an induced spin density 



given by the Kubo formula 
rL 



^ a, a' 

{a\s''{x,y)\a'){a'\jy\a) 
Ea - Ea' + if] 



)-f{Ea') 



En' 



(2) 



where s^{x^y) is the z-component of the local spin den- 
sity operator (we assume spin 3h/2 for the heavy- holes), 
jy = edH/dhky is the longitudinal electrical current op- 
erator for holes, and Ea and \a) are the eigenenergies 
and eigenstates of the full Hamiltonian H which includes 
the SO coupling term (^), the disorder potential, and 
the hard-wall confining potential, rj must be chosen to 
be small compared to all intensive energy scales, such 
as Ep^ h/r^ and the characteristic SO coupling strength 
Aso = 2Xkp, but larger than the finite system level spac- 
ing. A detailed discussion of subtleties related to the 
proper choice of r] can be found in Ref. 17]. 

Typical results of the calculated spatial profile of spin- 
density induced by a longitudinal electric field Ey are 
shown in Fig.^a). For comparison we show in Figs.^b) 
the transverse spin current, also obtained using linear 
response theory lEO^E^- both panels, Aso/Ef = 
0.4, {H/t)/ Ep = 0.3, and the diffusive regime is obtained 
by assuming system size L = 50 [/c^^] which is larger 
than the mean free path, / = 6[/c^^]. The z-component 
of the induced spin polarization has maxima near the 
strip edges, while the spin current has a maximum in the 
middle of the sample and decays toward the edges. We 
note that the peak value of the spin current divided by 
the applied electric field is comparable to the spin Hall 
conductivity calculated in a previous study with peri- 
odic boundary conditions rather than hard walls in the 
x-direction 17]. 

In Figs, ^c) and (d) we illustrate the disorder depen- 
dence of the peak value of the edge spin density, S^^^^^ 
and of the transverse spin current in the middle of the 
strip, j:! 

bulk' With increasing disorder strength, both 



edge 



and j 



x,bulk 



decrease gradually and vanish in the 
strong scattering limit, h/r ^ Aso- The striking simi- 
larity of the dependences on disorder strength of the two 
quantities led us to a numerical observation that. 



VF S, 



edge ^ Jbulk' 



(3) 



where vf is the hole Fermi velocity. Introducing the char- 
acteristic SO coupling length Lso = vpTrh/Aso and time 
tso = ^tt/Aso, we can rewrite Eq. as 



(4) 



Here we emphasize that both Lgo and tgo have an intrin- 
sic, disorder independent origin. Eq. suggests that 
the spin polarization accumulated in the region of size 
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edge ' 



is proportional to the amount of incoming 
spin current from the middle of the sample in a time scale 

^so • 



3 



^ a) SO dependence 










Aso/Ef=0.3 


3- y 


Aso/Ef=0.5 


R 


Aso/Ef=0.6 





X [kp-'; 



FIG. 2: Spatial profile of the 2;-component of accumu- 
lated spin density as a function of spin-orbit coupling for 
{h/T)/EF = 0.3 and L = AAkp^ (a), the disorder strength 
for Aso/Ef = 0.4 and L — 4:4kp^ (b), and of the system size 
for {n/r)/EF = 0.3 and Aso/Ef = 0.4 (c). System size de- 
pendence in a 2DEG for {n/T)/EF = 0.15 and Aso/Ef = 0.4 
(d). 



To further examine the validity of this ansatz we 
plot in Figs.[2Ia) and (b) the dependence of the spin po- 
larization profile on Ago and on h/r. As seen in Fig.I2fa), 
the spin density builds up more gradually from the cen- 
ter towards the edge of the strip for weaker SO coupling. 
Fig-E^b) demonstrates explicitly that the envelope of the 
spin polarization curve changes only weakly with disor- 
der. Calculations presented in Fig. [2|^c) further demon- 
strate that the edge spin accumulation properties do not 
change once the system size exceeds the SO length Lgo- 

Our findings on the disorder dependence of spin- 
accumulation differ from expectations based on analo- 
gies to other regimes. In the case of strong disorder the 
accumulation time scale is the spin-relaxation time 
Ts ~ Aj^r~^ while the accumulation length scale is the 
spin-diffusion length Ig = ^/Dts where D is the diffu- 
sion coefficient. For a given spin Hall current, analogies 
with this regime predict accumulation on the spin dif- 
fusion length Is 05 ^7 or the spin-precession length 
scale I/so [23, which corresponds with Ig in the dirty 
limit [5^, and predict that the magnitude is proportional 
to Ts oc . The applicability of this analogy in the 
strong SO scattering limit is not established theoreti- 
cally, however, and our numerical results suggest that 
it is doubtful. Consistency is not achieved by simply re- 
placing the weak SO interaction expression for the spin 
relaxation time by the strong SO interaction result that 
Ts ~ r, which would imply that the spin accumulation is 
proportional to r0. [T^l2n|. 

We can understand why the Kubo formula predicts 
that the induced spin-density is weakly dependent on 
scattering rate in the diffusive limit by considering the 
matrix elements of the spin-density and current opera- 



tors in Eq. Because the paramagnetic system is even 
under time-reversal and the spin-density is odd, it fol- 
lows that its static linear response originates from the 
odd dissipative term in the Kubo formula 



dy V- 9f{Ea) 



dEa 



where 6rj{E) 



{a\s'{x,y)\a'){a'\jy\a)S^{Eo. - Eo^.), (5) 

7r[£; 2+772] Typical properties of these ma- 
trix elements are illustrated in Fig. |31 where by "ME" we 
denote the above product of the spin and current matrix 
elements. As in all transport properties, current oper- 
ator matrix elements in the absence of disorder play a 
key role. For 2DHG strips, the component of momen- 
tum along the strips is a good quantum number whereas 
the transverse momentum is replaced by discretestand- 
ing wave or one-dimensinal subband indices. Both spin 
and charge operators have matrix elements that are di- 
agonal in the ID wavevector and in subband index. It 
follows that the induced spin-density that we calculate is 
(at T = 0) proportional to r]~^ in the absence of disor- 
der, diverging when the 77 ^ limit is taken. Unsurpris- 
ingly disorder, which limits the lifetime of these ID band 
states, must be considered to obtain a finite result. In 
the limit of the quasi ID limit or very weak disorder, the 
lifetime broadening by {h/r) of the band states is smaller 
than the subband separation, and it follows from Eq. (|5J 
that the induced spin-density will be proportional to r 
[23,123. This is the limit in which the ID subbands con- 
tribute independently to the induced spin density and 
can be treated as separate bands in a Boltzmann trans- 
port theory. 

The ID band limit applies only when h/r is smaller 
than the subband splitting, i.e., when h/r < hvp/L^ 
and therefore when the mean free path I = vpr \^ larger 
than As emphasized above, our calculations were 
performed in the regime in which / is small compared 
to the strip width. As illustrated in Fig. 13 intersub- 
band transitions contribute to the induced spin-density; 
electric field induced intersubband coherence cannot be 
neglected. Although the contribution from each transi- 
tion is proportional to r, the number of transitions that 
are incoherently mixed and contribute is proportional to 
r~^, leading to our finding of weak r dependence. 

We found that in the diffusive region, in which the 
system size is larger than the mean-free path, the matrix 
element products that appears in Eq. (0) at E^ — Ea' 
changes slowly with disorder. Exact disorder eigenstates 
at nearby energies mix ID band states with disorder free 
energies that vary over the range h/r and in the regime 
of interest strongly mix states belonging to different sub- 
bands. Eq. ^ shows that the induced spin-density is 
proportional to the quantity plotted in the main panel 
of Fig. OJa), which varies rapidly with scattering rate 
when the mean free path exceeds the system size, i.e. 
low values of {h/r)/EF^ and slowly at larger scattering 
rates as expected from the above argument. Fig. |3fb) 
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FIG. 3: Matrix elements that appear in Eq. © as a function 
of the scattering rate at Ea — E^' = (a) and as a function of 
the energy difference at {h/r)/EF = 0.14 (b). Aso/Ep = 0.3 
and L = 44A;^^ in both panels, (c) Shematic of the spin 
Hall accumulation in the strong SO coupling regime. Spatial 
profile of ^-component of spin density for (d) Aso/Ep = 0.6, 
0.4 (e), and 0.3 (f) at {h/r)/EF = 0.3. 



shows the eigenenergy diflPerence dependence of the ma- 
trix element product which remains large over the energy 
interval Aso- 

The above detailed analysis suggests the following 
physical picutre of the spin Hall accumulation in the 
strong SO interaction regime. Here we consider a re- 
flection process from a state with a momentum to 
another state k at the left boundary as shown in the 
Fig. 3. In the following we assume the width of stripe L 
is infinitely large, and focus on the major band. Note 
that the SO Hamiltonian has a form Hgo = —^0^) ' 
where the effective Zeeman field is given as A(k) = 
2A/c"(— sin(a(/)), cos(a(/)), 0) with a = 3 for 2DHGs and 
a = 1 for 2DEGs and = tdin~^ {ky/kx). Thus in a 
steady state the spins are parallel to the direction of this 
effective field A = A/|A|. The dynamics of the direc- 
tion of spins m = S/S is described by the equation of 
motion 

-m(k, t) = m(k, t) X — ^ + f __ . Vk j m(k, t) 

(6) 

in the clean limit. The spin precesses around A(k') until 
the reflection and changes its direction of the precession 
axis to A(k). In the presence of the electric field z- 
component spins are generated as discussed in the orig- 
inal paper of the intrinsic spin Hall effect 8]. Assum- 
ing m(k, t) A(k) + ^m(k, t), where Sm is linear in 
eE, we obtain {d/dt)'^Sm^{k,t) = -\A/h\'^{Sm^{k,t) - 
m^(k, E)) from Eq.©. Here note that the static compo- 
nent 

m^(k,E) = aeE • [z x k]/A:|A(k)| (7) 



causes the spin current density = 
^Ek,M^±/(^k,M)(^^?(k,E)^) = a,^[Exz] 

in the bulk 8], where s = 3/2 for 2DHGs. Note 
that the initial state with k' has a negative z- 
component of spin = —aeE\k'^\/k^\A{k.')\, while 
the right moving reflected state with k has a pos- 
itive = ae£^|/ca;|//c| A(k)|. A typical amount of 
z-component of spin for all left moving states {k'^ < 0) is 
^Efei<o,fe;,M/(^^k',M)^"^?(k'>E) ^ -saeE/AnvF = 
— Ij^I/'^f, and the same amount with opposite sign 
Mj^\/^F is for all right moving states {kx > 0). There- 
fore in the bulk of the sample these components cancel 
each other. On the other hand, near the edge this 
cancellation cannot be accomplished because it takes 
a finite time to change the averaged z-component spin 
from m^(k\ E) to m^(k, E) by the reflection as indicated 
schematically in Fig. 3(c). The net spin density arises 
from this edge induced coherence between different 
Bloch states and falls to zero on the length scale Lso'^/tt 
because of averaging over scattering angles as 

5,(x) = -^2expf-5-^). (8) 

Here the factor 2 /it comes from averaging of scattering 
angles. As shown in Fig.3(d)-(f) this ansatz is in a very 
good agreement with our numerical Kubo formula re- 
sults. Here we shifted the x coordinate as X ^ X «^edgG5 
where Xedge is the point of the peak of 5'^. These results 
show a qualitative correspondence between spin accumu- 
lation and the bulk SHE with the standard spin-current 
definition 0. 

We end this theoretical sections by noting that qualita- 
tively different behavior is expected for two-dimensional 
electron systems with Rashba SO coupling. Fig. I2fd) 
shows oscillations of the z-component of the 2DEG spin 
density whose amplitude decreases with increasing sys- 
tem size. The lack of clear edge spin accumulation 
in this system we obtained numerically is consistent 
with the suppressed bulk spin current of this model 
0, O Q El • A clear edge spin accumulation 
could appear in the non-linear transport regime with a 
strong electric field in a ballistic nanostructre.|2^ 



III. EXPERIMENT 

We now show that the above theoretical phenomenol- 
ogy of the SHE in strong SO coupling regime is consis- 
tent with our experimental data in modulation doped, 
GaAs/AlGaAs 2DHG devices. To demonstrate the inde- 
pendence of the SHE on the channel width we compare 
measurements in the 1.5 /im channel 4] and in a new de- 
vice with a 10 /im wide channel, lithographically defined 
in the same 2DHG system with / ^ 50 nm. In Fig. HJa) 
we show schematically the experimental setup with the 
LED current, I led, driving the electroluminescence at 
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FIG. 4: Schematic configuration of lateral p-n junction to 
detect spin accumulation (a): 2DHG channel bordered by 
an n-type region which forms LED. Driving a current I led 
through the LED results in electroluminescence indicated by 
the yellow line. Driving a current Ip along the 2DHG chan- 
nel induces the SHE spin accumulations at the edges of the 
2DHG. Light emission from the p-n junction recorded by a 
CCD camera (b) . Electron microscope image of the microde- 
vice with symmetrically placed LEDs at both edges of the 
2DHG channel (c). Emitted light polarization of recombined 
light in the p-n junction, (d) and (e), for the current flow 
indicated in (c). 



Fig. [pimply an approximate general form of the dimen- 
sionless edge polarization, 



edge 

hn 



9. 



(9) 



where gs is a numerical factor. In the measured system 
the 2DHG density n = 2 x lO^^cm-^ vf = 10^ms-\ 
^so/Ef ~ 0.4, and {h/T)/EF ^ 0.1. The system is 
therefore in the strong SO coupling regime and the cor- 
responding numerical factor gs = 7.4. From Eq. Q 
we obtain that the expected polarization for the mea- 
sured sample reaches ~ 8% and the width of the accu- 
mulation area is of order Lso ~ lOnm. The theoreti- 
cal ~ 8% polarization is consistent with the measured 
value of ~ 1% assuming an effective LED recombination 
width of the order ~ lOOnm. This number cannot be 
precisely determined experimentally because of the reso- 
lution limit set up by the wavelength of the emitted light 
(~ 800nm). However, comparisons between experiments 
in the 1.5 /im and 10 /im 2DHG channels, analysis of 
the digital images of the active p — n junction area, and 
simulations of device I — V characteristics *^ confirm a 
sub-micron width of the recombination region near the 
p — n junction. 



the edge of the hole strip, and the 2DHG channel cur- 
rent, /p, inducing the SHE edge spin accumulation. Dig- 
ital camera image of the light emitted from the p — n 
junction is shown in Fig. ^b). Fig. ^c) is a scanning 
electron microscope image of the whole microdevice with 
symmetrically placed LEDs at both edges of the 10 /im 
channel. Circular polarization of the spectral peak cor- 
responding to recombination of electrons with the spin- 
polarized 2DHG states at the edges of the channel are 
shown in Figs, ^d) and (e). For more details on our 
co-planar LED devices and on the spectroscopical analy- 
sis of the electroluminescence data see Refs. 0,113 and 
citations therein. 

Consistency of the measured data with the basic SHE 
phenomenology is demonstrated by the detected oppo- 
site sign of the z-component of spin accumulation at op- 
posite edges (compare red and blue curves in Fig. ^d) 
or (e)) and by the observed flipping of the sign of spin 
polarization at a given edge upon reversing Ip (com- 
pare, e.g., red curves in Figs. 2|Jd) and (e)). The mag- 
nitude of the Ip current in the 10 /im channel was ad- 
justed so that the corresponding longitudinal electric 
field, Ey ^ 0.15V/im~^, was comparable to the field used 
in the experiment in the 1.5 /im channel 3- Spin po- 
larizations of order ~ 1% observed in both devices con- 
firm the expected independence of the SHE signal on the 
width of the diffusive channel. 

We now compare experimental and calculated magni- 
tudes of the edge spin polarization. Numerical data in 



IV. SUMMARY 

In summary we have studied the edge spin accumu- 
lation due to the SHE in a 2DHG both theoretically 
and experimentally. We found theoretically a clear cor- 
respondence between the edge spins and the bulk spin 
currents. The calculated spin-density is nearly indepen- 
dent of the disorder scattering rate in magnitude and 
is localized within a spin-precession length of the edge. 
The magnitude increases with the strength of the SO cou- 
pling. For the parameters used in the experimental study, 
we predict an 8% spin-polarization at the edge. These 
predictions are consistent with the experimental finding 
of ~ 1% optical polarization averaged over a distance 
~ lOLso- The experiments also confirmed that opposite 
spin accumulations at the edges can be be separated over 
large distances in diffusive conduction channels and that 
their magnitude is not affected by the channel width. 
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